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Overlap Hypercube Fermions in QCD 
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We present simulation results obtained with overlap hypercube fermions in QCD near the chiral limit. We relate 
our results to chiral perturbation theory in both, the e-regime and in the p-regime. In particular we measured the 
pion decay constant by different methods, as well as the chiral condensate, light meson masses, the PCAC quark 
mass and the renormalisation constant Za- 



in 
o 



in 
o 



X 



1. INTRODUCTION 

We report on our recent simulation results 
for QCD on the lattice with chirally symmetric 
fermions. The aim is to relate our numerical re- 
sults to chiral perturbation theory (^PT), as an 
effective theory of the strong interaction at low 
energies. The simulation of QCD, as the funda- 
mental theory, allows for a determination of free 
parameters in this effective theory. 

In Section 2 we review the construction of the 
hypercube fermion (HF) and show recent appli- 
cations at finite temperature. We then proceed 
to the corresponding overlap hypercube operator, 
and we present results on its locality at strong 
gauge coupling. 

In Section 3 we turn our interest to the e- 
regime, which is the main focus of this work. In 
this setting, the linear size of the volume is shorter 
than the correlation length (which is given by the 
inverse pion mass). As an important virtue, eval- 
uations in this regime do not need an extrapola- 
tion to the infinite volume: physical quantities — 
in particular the Low Energy Constants (LEC) 
of the chiral Lagrangian — can be identified di- 
rectly in a small box, with their values in infinite 
volume. 

In Section 4 we add results obtained in the 
(standard) p-regime, where the xPT works by 
expanding in the momenta of the light mesons 
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involved, and in Section 5 we arrive at our con- 
clusions. 

Throughout this work we refer to quenched 
simulations using the (standard) Wilson gauge 
action, and we deal with a physical volume of 
V = (1.48 fm)^ X 2.96 fm at two lattice spacings 
a (measured by the Sommer scale pQ). Prelimi- 
nary results were presented before in Refs. |2|. 

2. OVERLAP HYPERCUBE FERMIONS 

For free fermions, perfect lattice actions are 
known analytically 3 . The corresponding lattice 
Dirac operator can be written as 



Dx,y = lt,Pt,{x -y) + \{x - y) , 



(1) 



with closed expression for and A in momentum 
space. Such actions emerge from iterated renor- 
malisation group transformations, hence they are 
free of any lattice artifacts regarding scaling and 
chirality. Dx,y is local (for non-vanishing A) but 
its range is infinite. Hence we tuned the renor- 
malisation group parameters for an optimal lo- 
cality of Dx,y] then we truncate it by means of 
periodic boundary conditions so that the sup- 
ports of Pfj,{x — y) and A(a; — y) are contained 
in {\xi, ~ yv\ < 1} for v = \ . . .A (in lattice units, 
a — 1). This truncation yields the free Hypercube 
Fermion (HF) which still has excellent scaling and 
chirality properties f5l. 

The HF is gauged by (normalised) sums over 
the shortest lattice paths between its sites. The 
paths correspond to products of fat link gauge 
variables, which are simply built as U^{x) — > 
(1 — a)U^{x) + ^ X] staples. Finally the links are 
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Figure 1. The spectral function crps, as a func- 
tion of the Matsubara frequencies uj, at = oo 
(on top, free case) and at finite Tc (below, inter- 
acting case) J^. These results are obtained from 
lattice data using the Maximum Entropy Method, 
as suggested in Ref. J^. 

amplified to restore criticality and to minimise 
the violation of the Ginsparg- Wilson relation. 

Due to the truncation and the imperfect gaug- 
ing procedure, the scaling and chirality are some- 
what distorted. Still the HF has highly favourable 
features, particularly in thermodynamics j4l5j . 
As an example, we show in Fig. ^ recent re- 
sults for the spectral function crps, which reveal a 
continuum-like behaviour up to much larger en- 
ergies than it emerges for the Wilson fermion. 

Exact chirality can be recovered by inserting 
the HF into the overlap formula 

= p{l + A/VA^) , A:=Do-p, (2) 



where p > 1 and Dq is some lattice Dirac operator, 
with Do = 75I?q75 (75-Hermiticity). 

• The standard formulation |S] inserts the Wil- 
son Dirac operator, Dq — I?W: which is then dras- 
tically changed in eq. ^ . We denote the resulting 
standard overlap operator, or Neuberger operator, 
as Z)ov— w- 

• Here we mainly study the case where instead 
the HF is inserted in the overlap formula 10), 
^0 — -Dhf El- This yields the operator Dqv-hf, 
which describes the overlap-HF. 

In both cases one arrives at exact solutions 
to the Ginsparg- Wilson relation jjnii and there- 
fore at an exact, lattice modified chiral symmetry 
im.^ However, in contrast to I?W: D^if is ap- 
proximately chiral already, hence its transforma- 
tion by the overlap formula Dhf -Dqv-hf, 
is only a modest modification. Therefore, the 
virtues of the HF are essentially inherited by the 
overlap-HF. Compared to Dov-w, the degrees of 
locality and of approximate rotation symmetry 
are highly improved |9ll4ll5in)j . If the gauge 
coupling is not too strong, the maximal impact 
of a unit source on a sink, as well as the maximal 
violation of rotation symmetry, decay exponen- 
tially in the distance. For instance, at (3 — 6 the 
exponent of these decays is almost doubled for 
the overlap-HF compared to Dov-w ^M- the 
latter case, locality was first studied in Ref. \n\). 

Here we are going to show mostly results at 
(3 = 5.85, which corresponds to a lattice spac- 
ing of a ~ 0.123 fm. We build the fat link with 
a = 0.52 and amplify all link variables by a fac- 
tor u ~ 1.28 (in the vector term we multiply 
another factor of 0.96). This is similar (but not 
exactly equivalent) to choosing the parameter p 
in Dov-w] for Dqv-hf we fix p = 1. These pa- 
rameters are chosen as a compromise between op- 
timising the locality and the condition number of 
A^A. Fig. 121 (on top) shows that the locality is 
significantly improved upon Dov-w (at p = 1.6) 

In Fig. El (below) we illustrate the locality at 
strong gauge coupling: at /3 = 5.7 (a ~ 0.17 fm) 



^The correctness of the axial anomaly in all topological 
sectors has been verified for the standard overlap operator 
in Ref. 12 , and for the overlap HF in Ref. 
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Figure 2. The locality of different overlap 
fermions, measured by the maximal impact of a 
unit source ipx on ipy at a distance r = ||a; — y||i. 
At the same value of f3, the overlap-HF is clearly 
more local, as the plot on top shows for (3 — 5.85. 
Below we see that its locality persist up to [3 — 5.6, 
where no exponential decay is visible for the stan- 
dard overlap formulation. 

the standard overlap operator is still barely local, 
if one chooses the optimal value p — 1.8. Simi- 
larly we optimise m = 1.6 for the HF and we still 
find a clear locality at this gauge coupling — this 
locality is in fact still stronger than the one ob- 
served for Dov-w a,t P — 6 (and p = 1.4, which 
is locality-optimal in that case). 
If we proceed to (3 ~ 5.6, the exponential decay 
evaporates for Dov-w, hence in this case the stan- 
dard overlap formulation does not provide a valid 
Dirac operator any more. On the other hand, if 
we insert the HF at m = 1.7 we still observe local- 
ity. Thus we see that the overlap-HF formulation 
provides chiral fermions on coarser lattices. 



3. EVALUATION OF THE LEADING 
LOW ENERGY CONSTANTS IN THE 
e-REGIME 

In the e-regime jl8) the correlation length ex- 
ceeds the box length, l/niT^ > L, and the obser- 
vables strongly depend on the topological sector 
|19| . Squeezing pions into such a tiny box may 
not be a realistic situation. Still, there is a strik- 
ing motivation for studying this regime: it allows 
for an evaluation of the LEG with their values 
in infinite volume, i.e. with their physical values. 
(Unfortunately, quenching brings in logarithmic 
finite size effects [20], but there is a window of 
volumes with decent values nevertheless). 

We first take a look at the overlap-HF indices: 
for 1013 configurations we obtained the histogram 
in Fig.|ni(on top). Fig.|Hl(below) shows our results 
for the topological susceptibility, compared to the 
continuum extrapolation of Ref. j21| . which used 
I?ov-w- Our susceptibilities are somewhat larger, 
but a rough continuum extrapolation agrees well 
with Ref. ED. 

3.1. Link to the chiral Random Matrix 
Theory 

Chiral Random Matrix Theory (RMT) conjec- 
tures the densities pn\z) of the lowest Dirac 
eigenvalues A in the e-regime [22], where z := 
XTiV . n = 1, 2, . . . numerates the lowest non-zero 
eigenvalues and v is the fermion index, which is 
identified with the topological charge ^I]. This 
RMT conjecture holds to a good precision for 
the lowest n and if L exceeds some lower 
limit (about 1.1 ... 1.5 fm, depending on the ex- 
act criterion) |2S]. Then the fit determines the 
scalar condensate E. From the cumulative den- 
sity of the first eigenvalue in the topologically 
neutral sector, the fit implies a chiral condensate 
of E = (299(8) MeV)3 for the overlap HF, in good 
agreement with the result of S = (300(7) MeV)^ 
for the standard overlap operator, see Fig. ^ 

3.2. Zero-mode contributions to the pseu- 
doscalar correlation 

We stay in the chiral limit and evaluate the 
time derivative of the zero-mode contribution to 
the correlation of the pseudo-scalar density P = 
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Figure 3. Histogram of overlap- HF indices (on 
top), and various results for the topological sus- 
ceptibility (below). 



with 



Di'!)\i;x,t) = 



(3) 



Tliis quantity was computed analytically to next- 
to-next-to-leading order (NNLO) in quenched 
xPT in Rcf. 24 , which also presented a numeri- 
cal study with D^v-w- We extend this study by 
including the overlap-HF. 

To NNLO the quenched chiral Lagrangian 
reads 



C 



NNLO 



-Str(AfJ7 + h.c.) 



'"-0 (^2 , "0 ^ .2 



-iK<i>oStr{MU - h.c. 
where U G SU{Nf\Nf)®SU{Nf\Nf)/SU{Nf\Nf), 




standard overlap 
overlap HF 
RMT prediction 



Figure 4. The cumulative density of the first 
Dirac eigenvalue Ai in the sector v = — stere- 
ographically projected and re-scaled with Y.V at 
optimal S — compared to the RMT prediction. 

"Str" is the corresponding "supertrace" , and Af 
is the diagonal matrix of fermion masses. The 
peculiarity of the quenched chiral Lagrangian is 



NNLO 



contains 



the presence of a scalar field $o- ^ 
five quenched LEG: Ft^, E, K, mo and ao. 

We consider the first term in the Taylor expan- 
sion of C|^|(i) around T/2. In a box x T it 
takes the form 



L2 



= + - (4) 



T 



T 



with s := t-T/2 and a := ao-2NcKF.^/Y, . The 
dependence on (i^^) originates from the Witten- 
Veneziano formula |25| . and mp was counted in 
Ois). 

We analysed the zero-modes that we collected 
at /5 = 5.85 on a 12^ x 24 lattice, and at /5 = 6 
on 16"^ X 32. In both cases, the physical volume 
amounts again to (1.48 fm)'^ x 2.96 fm. This is 
expected to be sufficiently large, referring to our 
experience with the chiral Random Matrix The- 
ory [231 and to the behaviour of the axial corre- 
lators in the e-regime (cf. Subsection 3.3). Our 
statistics is listed in Table We fixed the value 
of {v'^) from our simulations, since it can deviate 
significantly from the continuum extrapolation, 
as Fig. O (below) shows. 

Then we performed a combined two parameter 
fit of the leading term in eq. (^J to our data in the 
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10.8 


221 


192 


80± 14 


-17± 10 


standard ov. 12^ x 24, /3 = 5.85 


1.6 


11.5 


232 


180 


74± 10 


-19 ±8 


standard ov. x 32, (3 — 6 


1.6 


10.5 


115 


94 


75 ±24 


-21 ±15 



Table 1 

Statistics and results for F^^ and a, based on the zero-mode contributions to (PP) ■ The values for F^ 
and a in the last two columns correspond to the fitting range Smax = 1 • 



topological sectors — 1,2 over the intervals s £ 

[-Smax, Smax], with Smax 1 . . . 3. Our rCSultS 

are shown in Fig. [S] and in Table ^ Our values 
for and a are a little lower than the results 
reported in Ref. (for the Neuberger operator 
in isotropic boxes). 
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Figure 5. The Low Energy Constants _FV and a, 
determined from the zero-mode contributions to 
(PP) , against the width of the fitting range Smax . 

However, we see that our values of F^^ and a 
agree well for different overlap fermions and for 
two lattice spacings. This result for F.^ is close 



to the physical value — extrapolated to the chi- 
ral limit — of 86 MeV but below typical 

quenched values, see Subsection 3.3. 

3.3. Axial vector correlations 

In quenched xPT, the axial vector correlator 
depends to leading order only on the LEG S and 
F^ |27|. The_prediction for (^14(^)^4(0)) (with 
Aj^it) := ^/'(x, t)7574-!/'(af, t), i.e. the bare ax- 
ial current at p = 0) is a parabola with a mini- 
mum a.t t = r/2, where F^/T enters as an ad- 
ditive constant. In a previous study we observed 
that L should again be above 1 fm for this pre- 
diction to work, and that the history in the sector 
V = may be plagued by spikes, related to the oc- 
currence of very small non-zero eigenvalues |28| . 
(This leads to the requirement of a large statis- 
tics, but a method called Low Mode Averaging 
was then invented as an attempt to overcome this 
problem UHl-) In the non-trivial sectors, E can 
hardly be determined, but can be evaluated 
|2I30| . Fig. shows our current results with three 
(bare) quark masses niq in the e-regime, in the 
sector = 1 and 2. The (flavour degenerated) 
mass is incorporated in the overlap operator as 

^ov(m,) = (l-^)i^(";+m, . (5) 

Our result is based on 45 propagators for each 
mass in the sector = 1, and 10 propagators 
for each mass in the sector \v\ — 2. Inserting the 
RMT result of S ~ (299 MeV)^, as well as the 
renormalisation constant Za ~ 1.17 (see Section 
4), a global fit leads to = (113 ± 7) MeV, in 
agreement with other quenched results | |29l24l30| . 

4. RESULTS IN THE p-REGIME 

At last we present our results in the p- 
regime, which is characterised by a box length 
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Figure 6. Results for the axial vector correlator at 
arUq — 0.001, 0.003 and 0.005 in the sectors = 
1 (red, larger around the centre) and 2 (blue), and 
fits to the formulae of quenched xPT. ( Our quark 
masses correspond to T,Vmq — 0.83, 2.5 and 4.2.) 

L ^ l/^TT-TT, so that the p-expansion of xPT [HT] 
apphes. We consider once more /3 = 5.85, a lat- 
tice of size 12^ X 24 and the bare quark masses 
aruq = 0.01, 0.02, 0.04, 0.06, 0.08 and 0.1 (in 
physical units: m, = 16.1 MeV . . . 161 MeV). 
We computed 100 overlap-HF propagators, and 
we first evaluated m^r in three different ways: 

(1) From the pseudoscalar correlator (P(t)P(O)). 

(2) From the axial vector correlator (^14(^)^44(0)). 

(3) From (P(i)P(O) - S{t)S{0)), where S = ^ip. 
The subtraction of the scalar density is useful 
at small rrig to avoid the contamination by zero 
modes. 

The results are shown in Fig. |3 (on top): they 
follow well the expected behaviour arn^ cx ruq, in 
particular for {PP — SS) . Our smallest pion mass 
in this plot, ~ 289 MeV, has Compton wave 
length Ri L/2, hence this point is at the edge of 
the p-regime. 

Fig. (below) shows the corresponding results 
for the vector meson mass, with a (linear) chiral 
extrapolation to nip = 1017(39) MeV. Generally, 
quenched results tend to be above the physical 
p-mass. Our result is in agreement with a simi- 
lar study using Z?ov-w [121 1 but somewhat above 
the (quenched) world data for rup 133^ (which are 
mostly obtained with chiral symmetry breaking 
lattice fcrmions). 
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Figure 7. On top: The pion mass evaluated in 
various ways. Below: the p-meson mass. 

In Fig. |S1 (on top) we show the PCAC quark 
mass obtained from the axial Ward identity, 



mpcAC 



Es{dtAl{x)Pm 

2[E.-(^t(^mo)>] 



(6) 



We observe also here a behaviour which is nearly 
linear in rUq , with a chiral extrapolation to 
ampcAcimq = 0) = -0.00029(64). 
Remarkably, mpcAC is close to ruq, which is not 
the case for Dov-w fS^*. Consequently the renor- 
malisation constant Za — ruq /mpcAC is close to 
1; it has the chiral extrapolation Za — 1.17(2), 
which has been used already in Subsection 3.3. 
This is in contrast to the large Za factors ob- 
tained for Dov-w^ at (3 = 5.85 (p = 1.6), 
Za — 1.45 was measured |32I30| . and at /3 = 6 
(p = 1.4) even Za 1.55 E3|. 

Finally we reconsider the pion decay constant, 
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Figure 8. The PC AC quark mass of eq. and 
the renormalisation constant Za — 'niq/m-pcAC- 



now evaluated in the direct way, according to 
2m, 



\PW) 



(7) 



by using either PP or PP — SS, see Fig. 
El The hnear extrapolation to rrig — 
yields F^,pp = 111.5(2.5) MeV, resp. 
F^^PP_ss = 104(9) MeV, which is compatible 
with the result in Subsection 3.2. In particular 
the behaviour of the {PP — SS) result for F^^ at 
small niq calls for an evaluation at yet smaller 
quark masses, which was reported in Section 3. 

5. CONCLUSIONS 

The HF operator is approximately chiral and 
it has favourable properties in particular in ther- 
modynamics. It can be inserted into the overlap 
formula, which yields the overlap-HF. This is a 
formulation of chiral lattice fermions with better 
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Figure 9. The pion decay constant from a direct 
evaluation in the p-regime, using the overlap-HF. 

rotation symmetry and locality than the stan- 
dard formulation. Therefore it provides access to 
chiral fermions on coarser lattices. 

In the e-regime we evaluated — at tiny quark 
masses — the pion decay constant from the axial 
vector correlator. Directly in the chiral limit we 
obtained results for the chiral condensate from 
the Dirac spectrum. Moreover, from the zero- 
mode contribution to the pseudoscalar correlator, 
we found again values of F^^ and the parameter 
a (a quenching specific LEG). The results are 
very similar for two lattice spacings and different 
overlap operators. Our result for Ft^ obtained 
with this method (in the given volume) is close 
to the physical value (in the chiral limit), but 
below other quenched values. 
We add that a topology conserving gauge action 
could be helpful in that regime |SS|- 

In the p-regime we obtained results for m^, 
mp and F^^. Compared to the standard overlap 
fermion, mpcAC is closer to the bare quark mass 
mq, hence Za is much closer to 1, which is pleas- 
ant for a connection to perturbation theory. Here 
we measured F^^ directly, and the chiral extrap- 
olation of this result agrees (within the errors) 
with the value obtained in the e-regime from the 
axial vector correlator. 



•^The two methods applied for the axial vector and the 
pseudoscalar correlator differ by a subtlety in the quenched 
counting rules for the e-expansion. 
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